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This letter reports results from the MINOS experiment based on its initial exposure to neutrinos 
from the Fermilab NuMI beam. The rates and energy spectra of charged current f M interactions 
are compared in two detectors located along the beam axis at distances of 1 km and 735 km. With 
1.27 x 10 20 120 GeV protons incident on the NuMI target, 215 events with energies below 30GeV 
are observed at the Far Detector, compared to an expectation of 336 ± 14.4 events. The data 
are consistent with disappearance via oscillations with |Am3 2 | =2.74j^26 x 10~ 3 eV 2 /c 4 and 
sin 2 (26»23) > 0.87 (68% C.L.). 



PACS numbers: 14.60.Lm, 14.60.Pq, 29.27.-a, 29.30.-h 

There is now substantial evidence [JEHSIEEHH 
that the proper description of neutrinos involves a ro- 
tation between mass and flavor eigenstates governed by 
the 3x3 PMNS matrix [1 [To| . The parameters of this 
mixing matrix, three angles and a phase, as well as the 
mass differences between the three mass eigenstates must 
be determined experimentally. The Main Injector Neu- 
trino Oscillation Search (MINOS) experiment has been 
designed to study the flavor composition of a beam of 
muon neutrinos as it travels between the Near Detec- 
tor (ND) at Fermi National Accelerator Laboratory at 
lkm from the target, and the Far Detector (FD) in the 
Soudan iron mine in Minnesota at 735 km from the tar- 
get. From the comparison of the reconstructed neutrino 
energy spectra at the near and far locations the oscilla- 
tion parameters |Am§ 2 and sin 2 (2# 2 3) are extracted. 

The Neutrinos at the Main Injector (NuMI) neutrino 
beam is produced using 120 GeV protons from the Main 
Injector. The protons are delivered in 10 /zs spills with 
up to 3.0xl0 13 protons per spill. The extracted protons 
are bent downward by 3.3° to point at the MINOS de- 
tectors. The global positioning system (GPS) defined the 
survey beam direction to within 12 m of the FD ^l| . Pos- 



itively charged particles produced by the proton beam in 
the 95.4cm long target (mainly 7r + and K + ) are focused 
by two pulsed parabolic horns spaced 10 m apart and are 
then allowed to d ecay in a 675 m long, 2 m diameter, evac- 
uated decay pipe [l2j . The proton beam [13J and tertiary 
muon beam [14| are monitored on a pulse-by-pulse basis. 
The target position relative to the first horn and the horn 
current are variable [15|. For the majority of the running 
period described here, the target was inserted 50.4 cm 
into the first horn to maximize neutrino production in 
the 1-3 GeV energy range. A total of 1.27xl0 20 protons 
on target (POT) were taken in this position and used for 
the oscillation analysis. The charged current (CC) neu- 
trino event yields at the ND are predicted to be 92.9% 
v^, 5.8% 77 M , 1.2% v e and 0.1% v e . The data described 
here were recorded between May 2005 and February 2006. 
The average livetime of the FD was 99.0% during this pe- 
riod. About one third of the total ND events provided a 
sufficiently large dataset for this analysis of ~10 6 events 
which were sampled throughout the run period. 

Both MINOS detectors ^| are steel-scintillator track- 
ing calorimeters 0] with toroidal magnetic fields aver- 
aging 1.3 T0. The steel plates are 2.54 cm thick. The 
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scintillator planes are comprised of 4.1 cm wide and 1 cm 
thick plastic strips. Each plane is oriented at 45° from 
vertical and at 90° with respect to its neighbors. The 
light from the scintillator strips is transported to the 
multi-anode photomultiplier tubes (PMT) by embedded 
1.2 mm diameter wavelength shifting (WLS) fibers. In 
order to cancel the majority of the uncertainties in the 
modeling of neutrino interactions and detector response, 
the two MINOS detectors are as similar as possible. For 
example, both detectors yield 6 — 7 photoelectrons (PEs) 
per plane for normally incident minimum ionizing parti- 
cles. However, the data rate in the ND is ~ 10 5 times 
larger than in the FD which has dictated certain design 
differences between them. 

The 5.4kton FD, 705 m underground, has 484 octago- 
nal, 8 m wide instrumented planes read out at both ends 
via Hamamatsu M16 PMTs 0. Eight WLS fibers from 
strips in the same plane, separated from each other by 
about 1 m, are coupled to each pixel. The coupling pat- 
tern is different at the two ends to allow resolution of 
ambiguities. 

The 0.98 kton ND, 103 m underground, has 282 irregu- 
lar 4x6 m 2 octagonal planes. Its geometry optimizes the 
containment of hadronic showers and provides sufficient 
flux return to achieve a magnetic field similar to the FD. 
Each strip is coupled via a WLS fiber to one pixel of a 
Hamamatsu M64 PMT 20]. The ND readout continu- 
ously integrates the PMT charges with a sampling rate 
of 53.1MHz to allow discrimination between successive 
Main Injector RF buckets. 

The data aquisition HH H3, accepts data above a 
threshold of 0.25 PEs. In the FD, the online trigger con- 
ditions require a hit within 100 (is centered on the time 
of the expected beam spill, at least 20 PEs inside a four 
plane window, or 4 hits in 5 consecutive planes. In the 
ND, all the data taken during the beam spill are retained. 
The trigger efficiency for both detectors exceeds 99.5% 
for neutrino events with visible energy above 0.5 GeV. 

The detectors are calibrated using an in-situ light in- 
jection system [24| and cosmic ray muons. LED gener- 
ated light signals are distributed to all the WLS fibers to 
track gain changes in the PMTs and electronics. The en- 
ergy deposited by through-going muons is used to equal- 
ize the response of all the scintillator strips. After cal- 
ibration, remaining time and position dependent varia- 
tions in the responses of the detectors result in an uncer- 
tainty in the relative energy scale between the two detec- 
tors of 2%. The overall energy scale for single hadrons 
and electrons was determined from the results of a test- 
beam experiment using a small, unmagnetized copy of 
the MINOS calorimeters (CalDet) j^. Stopping muons 
are then used to relate the results from CalDet to the re- 
sponse of the ND and FD. From these studies, the uncer- 
tainty on the absolute hadronic energy scale is estimated 
to be 6%. 

The simulation of the production and detection of neu- 



trinos commences with a model of hadron production 
in the target using FLUKA05 j2|| , which has uncertain- 
ties at the 20-30% level stemming from a lack of rele- 
vant thick target hadron production data. The shower 
products are transported through the horn focusing sys- 
tem and decayed in a GEANT3 27] simulation that in- 
cludes the horns, beamline material and the decay pipe. 
The neutrino event generator, NEUGEN3 |23|. is tuned 
to existing CC cross-section data where present uncer- 
tainties below 10 GeV are at the 20% level. The prod- 
ucts of the neutrino interaction are propagated out of 
the iron nucleus using the INTRANUKE |29( code from 
within NEUGEN3. Some of the energy of absorbed pi- 
ons is transferred to clusters of nucleons as motivated by 
Ref. The response of the detector is simulated us- 
ing GEANT3 with the GCALOR 31] model of hadronic 
interactions. The final step in the simulation chain in- 
volves photon generation, propagation and transmission 
through the WLS fiber and conversion to photoelectrons 
in the PMTs. 

In CalDet, GEANT3 with GCALOR is found to repro- 
duce the hadronic and electromagnetic (EM) responses of 
the detector to single particles to 4% and 2%, respec- 
tively. Below 10 GeV, the hadronic energy resolution 
was measured to be 56%/^E[GeV] © 2% [3^] and the 
EM resolution was measured to be 21A%/^E[GeV] © 
i.l%/E[GeV] H3. The muon energy resolution AE^/E^ 
varies smoothly from 6% for E M above 1 GeV where most 
tracks are contained and measured by range, to 13% at 
high energies, where the curvature measurement is pri- 
marily used. 

The initial step in the reconstruction of the FD data 
is the removal of the eightfold hit-to-strip ambiguity us- 
ing information from both strip ends. In the ND, tim- 
ing and spatial information is first used to separate in- 
dividual neutrino interactions from the same spill. Sub- 
sequently, tracks are found and fitted, and showers are 
reconstructed, in the same way in both detectors. For 
CC events, the total reconstructed event energy is 
obtained by summing the muon energy and the visible 
energy of the hadronic system. 

The FD data set was left blind until the selection pro- 
cedure had been defined and the prediction of the unoscil- 
lated spectrum was understood. The blinding procedure 
hid a substantial fraction of the FD events with the pre- 
cise fraction and energy spectrum of the hidden sample 
unknown. Events are pre-selected in both detectors, by 
requiring total reconstructed energy below 30 GeV and 
a negatively charged track to suppress events that orig- 
inate from tt~ or . The track vertex must be within 
a fiducial volume such that cosmic rays are rejected and 
the hadronic energy of the event is contained within the 
volume of the detector. The event time must fall within 
a 50 /is window around the spill time. Cosmic ray back- 
ground is suppressed further in the FD by requiring the 
track to point within 53° of the neutrino beam direc- 
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FIG. 1: Data and tuned MC predictions for the PID variable 
in the ND (top) and FD (bottom). The arrows depict the 
positions of the ND and FD selection cuts. The FD MC dis- 
tribution for CC events uses the best fit parameters discussed 
in the text. 
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FIG. 2: Energy spectra in the MINOS ND for three of the 
six beam configurations before and after the 15 parameter 
beam tuning procedure. The target location was modified to 
produce the different spectra: a) nominal, b) target at 90 cm 
from nominal, c) target at 240 cm from nominal. The lower 
inset shows the ratio of data to MC before and after tuning. 



tion. The prc-sclcctcd event sample is predominantly 
CC with a 8.6% NC background estimated from Monte 
Carlo (MC) simulations. The fiducial mass of the FD 
(ND) is 72.9% (4.5%) of the total detector mass. 

A particle identification parameter (PID) incorporat- 
ing probability density functions for the event length, the 
fraction of energy contained in the track and the average 
track pulse height per plane provides separation of CC 
and NC events. The PID is shown in Fig.njfor ND and 
FD data overlaid with simulations of NC and CC events 
after the beam reweighting procedure described below. 
Events with PID above -0.2 (FD) and -0.1 (ND) are se- 
lected as being predominantly CC in origin. These values 
were optimized for both detectors such that the resulting 
purity of each sample is about 98%. The efficiencies for 
selecting CC events in the fiducial volume with en- 
ergy below 30 GeV are 74% (FD) and 67% (ND). From 
the absence of any events less than 20 [is before and less 
than 30 fis after the spill time, the remaining non-beam 
related background in the FD is estimated to be less than 
0.5 events (68% C.L.). Background from interactions 
in the rock surrounding the FD is estimated from MC 
to be below 0.4 (68% C.L.) events. The corresponding 
backgrounds in the ND are negligible. 

To constrain hadron production, a series of six runs 
of similar exposure was taken where the position of the 
target and the magnitude of the horn magnetic field were 
varied. Comparisons of the ND energy spectra with MC 
simulations, shown in Fig. [3 showed an energy dependent 
discrepancy that changed with the beam settings. This 
implied beam modeling, rather than detector or cross- 
section effects, was the primary cause. To bring the MC 
into better agreement with the data, a tuning of the beam 
MC was performed in which pion production off the tar- 
get was smoothly varied in transverse and longitudinal 
momentum with respect to the FLUKA05 input, as was 
the overall kaon yield. In addition, the potential sys- 



tematic effects of the beam focusing, NC background, 
Vp energy scale and offset were allowed to vary. All of 
these parameters were found to lie within two standard 
deviations of their nominal values. Fig.[2J shows the ef- 
fect of the full beam parameter tuning for the spectra 
corresponding to three different target positions. The 
resulting agreement is improved in all beams across the 
1-30 GeV neutrino energy region. 

The measurement of the energy spectrum at the ND 
is used to predict the unoscillated spectrum at the FD. 
The oscillation hypotheses are then tested relative to this 
prediction. The prediction must take into account the 
ND and FD spectral differences that are present, even in 
the absence of oscillations, due to pion decay kinematics 
and beamline geometry. These introduce a ND /FD shape 
difference of up to ~20% on either side of the peak. 

There are two distinct approaches to the beam extrap- 
olation. The ND Fit method focuses on minimizing the 
remaining ND data and MC differences by modifying MC 
parameters associated with neutrino interactions and de- 
tector response. The FD MC is then re-weighted with 
the best-fit values of these parameters. For the results 
presented in this paper, the Beam Matrix method [34| is 
used, in which agreement between MC and data is much 
less important because the ND data are used to measure 
all the effects common to both detectors, such as beam 
modeling, neutrino interactions and detector response. 
It utilizes the beam simulation to derive a transfer ma- 
trix that relates s in the two detectors via their parent 
hadrons. The matrix element Mij gives the relative prob- 
ability that the distribution of secondary hadrons which 
produce s of energy Ei in the ND will give s of 
energy Ej in the FD. The ND reconstructed event en- 
ergy spectrum is translated into a flux by first correcting 
for the simulated ND acceptance and then dividing by 
the calculated cross-sections for each energy bin. This 
flux is multiplied by the matrix to yield the predicted, 
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FIG. 3: Comparison of the Far Detector spectrum with pre- 
dictions for no oscillations for both analysis methods and for 
oscillations with the best-fit parameters from the Beam Ma- 
trix extrapolation method. The estimated NC background 
is also shown. The last energy bin contains events between 
18-30 GeV. 



unoscillated FD flux. After the inverse correction for 
cross-section and FD acceptance, the predicted FD visi- 
ble energy spectrum is obtained. 

In total, 215 events are observed below 30 GeV com- 
pared to the unoscillated expectation of 336.0±14.4. The 
error is due to the systematic uncertainties described be- 
low. In the region below 10 GeV, 122 events are observed 
compared to the expectation of 238.7±10.7. The ob- 
served energy spectrum is shown along with the predicted 
spectra for both extrapolation methods in Fig.EI 

Under the assumption that the observed deficit is due 
to v„ — > v T oscillations jH, HE , a fit is performed to 
and sin 2 (26*23) using the expres- 



the parameters |Am 2 2 | 
sion for the Vn survival probability: 



1 



sin 2 ( 26*23 ) sin 2 ( 1 . 2 7 A m\ 2 — ) 
E 



(1) 



where L[km] is the distance from the target, i?[GeV] 
is the neutrino energy, and |Arn 2 2 | 38] is measured in 
eV 2 /c 4 . The FD data are binned in reconstructed event 
energy and the observed number of events in each bin 
is compared to the expected number of events for this 
oscillation hypothesis. The best fit parameters are those 
which minimize x 2 =— 2 In A where A is the likelihood ra- 
tio: 
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where 0$ and ei are the observed and expected numbers 
of events in bin i, and the As|/cr 2 . are the penalty terms 
for nuisance parameters associated with the systematic 
uncertainties. The expected number of events depends 
on |Ato 2 2 |, sin 2 (26*23) an d the Sj. The choice of these 
systematic effects and their estimated uncertainties are 
described below. The e» include the small contribution 



from selected v T events produced in the oscillation pro- 
cess. 

The effects of different systematic uncertainties were 
evaluated by modifying the MC and performing a fit on 
this in place of the data. The differences between the 
fitted values obtained with the modified and unmodi- 
fied MC are listed in Tabled The largest effects are: 
(a) The uncertainty in the fiducial mass in both detec- 
tors, uncertainty in the event selection efficiency and the 
POT counting accuracy gives a 4% uncertainty on the 
predicted FD event rate, (b) The absolute hadronic en- 
ergy scale from a combination of test beam measurements 
and calibration accuracy is known to 6% as discussed 
above. This is added in quadrature to the uncertainty 
in the effect of intra-nuclear re-scattering estimated at 
±10% of the hadronic energy. The total hadronic energy 
scale uncertainty is therefore ±11%. (c) The NC com- 
ponent was varied in a fit to the PID data distribution 
in six energy bins in the ND. A 50% uncertainty was 
estimated to encompass the differences between the fit 
and NC MC. At the current level of statistics, uncertain- 
ties from CC cross-sections, muon momentum, relative 
ND/FD energy calibration, remaining beam uncertain- 
ties and reconstruction were found to be negligible. As 
an example, in the absence of any beam tuning, the best 
fit value only shifts by 0.2x 10" 5 eV 2 /c 4 . 

In fitting the data to Eqn.Q sin 2 (26*23) was constrained 
to lie in the physical region and the main systematic un- 
certainties ((a), (b) and (c) in TableHJ) were included as 
the nuisance parameters. The resulting 68% and 90% 
confidence intervals arc shown in Fig.H] as determined 
from A% 2 =2.3 and 4.6, respectively 39 ] . The best fit 
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is lAr 
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= (2.74^) x io-3 eV 2 / c 4 



and sin 2 (26 , 23 )>0.87 at 68% C.L. |40j with a fit prob- 
ability of 8.9%. At 90% C.L. (2.31<|Am§ 2 |< 3.43) x 
10" 3 eV 2 /c 4 , and sin 2 (20 23 ) > 0.78. The data and best 
fit MC are shown in Fig-EI At the best fit value, the 
MC predicts 0.76 v r events in the final sample. If the 
fit is not constrained to be within the physical region, 
\Am 2 32 1=2.72 xl0~ 3 eV 2 /c 4 and sin 2 (26»23)= 1-01, with a 
0.2 decrease in x 2 - With additional data, it is expected 
that the systematic uncertainties will be reduced. 
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(a) Normalization (± 4%) 


0.05 
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(b) Abs. hadronic E scale (± 11%) 
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(c) NC contamination (± 50%) 


0.09 


0.050 
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